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^P' ABSTRACT 

P ; 

We examine string (vortex) formation at a quench for a weakly-coupled global U(l) 
theory when the excitation spectrum is non-relativistic. It is so similar to vortex 
production in the corresponding relativistic plasma as to reinforce arguments for 
O^J ' the similarity of vortex production in the early universe and in low-temperature 

many-body physics. 

I „, 

structure of the universe can be attributed to cosmic strings formed at phase tran- 
sitions in the Grand Unification era. Unfortunately, given the unlikely event of 
observing a cosmic string (vortex) directly it is difficult to make the case com- 



C3 



1. Introduction 



§^ pelling. 

■ However, the formation of topological defects like vortices during symmetry- 

breaking phase transitions is not unique to the early universe but generic to many 
physical systems. In particular, recent experiments on the production of vortices 
in superfluid A He and 3 He a have excited considerable interest. A priori, we 
might expect the early universe to have little in common with many-body systems. 
Its highly relativistic regime is characterised by initial temperatures (or energies) 
T ^> m for all particle masses m (in units ks = c = 1), whereas low temperature 
non-relativistic media are characterised by T < m, to freeze out antiparticles, but 
with T ~ /x nr , the non-relativistic chemical potential. Nonetheless, the proposi- 
tion has been made by the experimental groups concerned that these, and other 
experiments!, may provide an insight on the early universe. In this they have been 
championed by some theoretical astroparticle physicists, most notably Zurek@, one 
of whose most recent articles0 was, indeed, titled Cosmological Experiments in Su- 
perfluids and Superconductors. 

In this talk I shall indicate how comparisons between vortex production in these 
very different environments can be more than superficial analogies by sketching out 
some preliminary quantitative steps. The talk has three components. 
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Firstly, I shall remind you of the generalities of vortex production and of the way 
a network of simple Z-vortices is characterised by field configuration probabilities. 
The next step is to show how, by varying the chemical potential, it is possible to 
interpolate between a relativistic scalar theory (which, in curved space-time, could 
be of the early universe), and a non-relativistic field theory of condensed bosons. 
Finally, these ideas will be given a concrete realisation in a model of Z-vortex 
formation by unstable long-wavelength Gaussian fluctuations. 

The reader seeking greater detail will find itin recent work by myself and my 
collaborators, Tim Evans and Alasdair Gillii'0'0, from which this talk is drawn. 

2. Vortex Production and Vortex Distributions 

What encourages us in the hope that vortex production in the early universe and 
the laboratory has close parallels is that the most plausible production mechanism 
refers to neither. We recapitulate it now for the simplest theory permitting vortices, 
that of a complex scalar field 0(x, t). The complex order parameter of the theory 
is (0) = r]e ia and the theory possesses a global 0(2) symmetry that we take to be 
spontaneously broken at a phase transition. The field <fi could be either a complex 
non-relativistic order field appropriate to a superfluid or a relativistic field in the 
early universe. 

Initially, we take the system to be in the symmetry-unbroken (disordered) phase. 
We have no reason to choose any particular initial field configuration, beyond the 
requirement that the field is distributed about = with zero mean. The simplest 
assumption is that, beginning at some time t = to, the 0(2) symmetry of the 
ground-state (vacuum) is broken by a rapid change in the environment inducing an 
explicit time-dependence in the field parameters. Once this quench is completed 
the 0-field potential V(<p) = — a\<f)\ 2 + &|0| 4 is taken to have the symmetry-broken 
form a > 0, b > of the familiar 'wine-bottle' bottom. The ground-state manifold 
(the circle S 1 , labelled by the phase a of (0)) is infinitely connected and the theory 
possesses global strings or vortices, labelled by a winding number n G Z. 

These strings cost considerable energy to produce. That they should appear at 
all follows from a general argument, due to Kibbled, which goes as follows. During 
the transition, the complex scalar field begins to fall from the false ground-state into 
the true ground-state, choosing a point on the ground-state manifold at each point 
in space. For continuous transitions, for which this collapse to the true ground-state 
occurs by phase separation, the resulting field configuration is expected to be one of 
domains within each of which the scalar field has relaxed to a constant ground-state 
value. 

If this is so, then the requirements of continuity and single valuedness will some- 
times force the field to remain in the false ground-state between some of the do- 
mains. For example, the phase of the field may change by an integer multiple of 2n 
on going round a loop in space. This requires at least one zero of the field within 
the loop, each of which has topological stability and characterises a vortex passing 



through the loop. The density of strings is then closely linked to the number of 
effective domains and the evolution of this density is, correspondingly, linked to the 
nature of the domain growth. When the phase transition is complete and there is 
no longer sufficient thermal energy available for the field to fluctuate into the false 
ground-state, the topological defects are frozen into the field. 

Consider an ensemble of systems evolving from one of a set of disordered states 
whose relative probabilities are known, to an ordered state as indicated above. If 
the phase change begins at time t then, for t > t , it is possible in principle to 
calculate the probability pt[$] that the complex field 0(x, t) takes the value 3>(x) at 
time t. Throughout, it will be convenient to decompose $ into real and imaginary 
parts as $ = -4j(<3?i +i<3? 2 ) (and <fi accordingly). This is because we wish to track the 
field as it falls from the unstable ground-state hump at the centre of the potential 
to the ground-state manifold in Cartesian field space. 

The calculation of pt[<&] will be performed later in our simple model. For the 
moment, consider it given. The question is, how can we infer the string densities 
and the density correlations from pt[3>]? That we can calculate them at all is a 
consequence of the fact, noted earlier, that the string core is a line of field zeroes. 
This is equally true for both relativistic and non-relativistic 0(2) theories. The 
zeroes of $ a (a=l,2) which define the vortex positions form either closed loops or 
'infinite' string i.e. string that does not intersect itself. Following HalperirJl3 we 
define the topological line density p(r) by 

p(r) = Y:Jds^[r-R n (s)]. (I) 

In (2.1) ds is the incremental length along the line of zeroes R n (s) (n=l,2,.. .) and 
is a unit vector pointing in the direction which corresponds to positive wind- 
ing number. Only winding numbers n = ±1 are considered, with higher winding 
numbers understood as describing multiple zeroes. If dAj is an incremental two- 
dimensional surface containing the point r, whose normal is in the jth direction, 
then Pj(r) is the net density of strings (i.e the density of strings minus the density 
of antistrings on dAj). 

Ensemble averaging at time t is understood as averaging over the field 

probabilities pt[&) as 

(F[$]) t = J V$p t [$] (2) 

In general this ensemble averaging is not thermal averaging since, out of equilibrium, 
we have no Boltzmann distribution. We shall only consider situations in which 

(Pi(r))t = 0, (3) 

i.e. an equal likelihood of a string or an antistring passing through an infinitesimal 
area. However, the line density correlation functions 



Cij(r;t) = ( Pi (r) Pj (0)) t 



(4) 



will be non-zero, and give information on the persistence length of strings. 
It follows that, in terms of the zeroes of 3>(r), Pi{r) can be written as 

Pi (r) = S^r^d^^d^r), (5) 

where 5 2 [$(r)] = 5[$i(r)]5[$ 2 (r)]. The coefficient of the 5-function in Eq.|] is the 
Jacobian of the transformation from line zeroes to field zeroes. It permits us to 
define a further line density that we shall also find useful, the total line density p(r) 

Pi(r) = S^^We^dj^dk^l (6) 

Unlike the case for p«(r) 

n(t) = (A(r)) t > (7) 

and measures the total string density in the direction i, without regard to string 
orientation. The isotropy of the initial state guarantees that n(t) is independent of 
the direction i. We note that the Jacobian factor multiplying the field (^-functions 
in Eq.[5] and Eq.[| guarantees that random field zeroes with no vorticity will not be 
counted. 

In general, the best that we can do is write (pi(r)) t as 

(p t (r)) t = J V$ p t [^ S'Mr^e^d^^dMr^ 

= J Va (leiM^ito^WI e ifd* aa * a)t (g) 

The simple model that we shall propose later assumes that p t [<&(r)] is Gaussian. 
In this case the system is solvable, as we shall now show. Let us assume that 

(*.(*)>* = o, 
($ a (r)9 i $ 6 (r)) t = 0, (9) 

and, further, that 

($ a (r)$ 6 (r'))t = W ab (\r-r'\;t) 

= 6 ab W(\r-r'\;t), (10) 

is diagonal. All other connected correlation functions are taken to be zero. Then 
all ensemble averages are given in terms of W(r;t), where r = |r|. In particular, 
(p~i(r))t separates as 

(*(r)) t = (5 2 Mv)}) t ({e^d^d^Di (11) 
or, equivalently, from Eq.|^ 

(p, t (r)) t = jva {\e ijk d^{r)d k ^ 2 {r)\) t (e*J ~ d ™°*°) t . (12) 



It follows!!, on first performing the a integration in the second factor, that 



n{t) 



1 

2^ 



W"(0;t) 



W(0;t) 

where primes on W denote differentiation with respect to r. Thus 

1 



(13) 



where £ is the length at time t that sets the scale in W(r; t). 
The density-density correlation functions 



(14) 



Cij(r;t) = A(r;t)6 i3 +B(r-t) 



8 



(15) 



have a much more complicated realisation^ in terms of W and its derivatives. We 
shall be looking for signs of anticorrelation, which enables us to determine the 
persistence length of strings in the network (i.e. how bendy they are). The bendier, 
the more string that will occur in small loops. This is important, since large scale 
structure formation in the universe requires a certain amount of infinite string. 

Nothing that we have said so far discriminates between vortices in a relativistic 
or a non-relativistic medium. This distinction appears in the definition of (...)*, to 
which we now turn. 



3. The Relationship Between Relativistic and Non-Relativistic Media 

In order to compare vortex formation in relativistic and non-relativistic media it 
is convenient to develop the formalism of the relativistic field to accomodate both. 
For simplicity we assume in each case that the system is initially in equilibrium in the 
disordered phase. The interpolation between relativistic and non-relativistic regimes 
is then effected by introducing the chemical potential fi into the relativistic theory, 
coupled to the conserved charge Q (particle number minus antiparticle number) 
arising from the 0(2) symmetry. Provided \x is small in comparison to a particle 
mass the introduction of such a potential will have little effect on a phase transition 
for the relativistic theory initiated by quenching the system from a high temperature 
T to, effectively, zero temperature. However, on increasing fi prior to the transition 
it becomes more costly to produce antiparticles and, if the initial temperature T is 
decreased to a value much less than /i, such antiparticles as are produced are frozen 
out. The system is then one of non-relativistic particles at a temperature much 
less than the particle rest masses. In this nonrelativistic regime the transition is 
induced by a quench in \i itself or equivalently, since \i determines the density, by a 
density (pressure) quench. 

There are two separate problems. The first is to show explicitly how the non- 
relativistic limit occurs. The second is to relate this non-relativistic limit (of a 



relativistic field) to a conventional a priori non-relativistic field theory of many- 
body physics. However, before we can do either we need to discuss non-equilibrium 
field dynamics. 



3.1. Field Dynamics with a Chemical Potential 

From the viewpoint given above the chemical potential is seen as determining the 
initial conditions for the subsequent dynamic^ fhrw hich we adopt the closed time 
path method (Schwinger-Keldysh formalism)t3'EjJi3li3. We have already assumed 
that, at the initial time to, we are in a disordered state with (0) = 0. Our ignorance is 
parametrised by the probability distribution p tQ [$] that, at time t , 0(t o ,x) = $(x). 
For the moment we take it as given. Whether we are in a relativistic or non- 
relativistic regime is largely encoded in Pt [$]. The subsequent, essentially generic, 
non-equilibrium field evolution is driven by a rapid change in the environment. 
Specifically, for t > to the action for the field is taken to be 

SM = J d'x^d^d^ - \m\t)<? a - \\{t){^. (16) 

where m(t), A(i) describe the evolution of the parameters of the theory under ex- 
ternal influences, to which the field responds. There is no penalty in using the 
Lorentz covariant Eq.|l6], even for non-relativistic media. If, initially, there are no 
antiparticles at low temperature, there will remain no antiparticles. As with $, it 
is convenient to decompose in terms of two massive real scalar fields O , a — 1, 2 
as — (0i +*02)/\/2, in terms of which S[(f>) shows a global 0(2) invariance, broken 
by the mass term if m 2 (t) is negative. 

The change of phase that begins at time to will, by the mechanism indicated 
earlier, lead to the appearence of vortices. We saw in the previous section that 
the vortex distributions at later times tf > t can be read off from the probability 
p t/ [$j] that the measurement of will give the value The evolution of Pt[&] 
from to to tf is most simply written as a closed time-path integral in which the field 
is integrated along the closed path C + © C_ of Fig.l, where = 0+ on C + and 
= 0_ on C-. If V(f)± = lla=i^0±,a and spatial labels are suppressed then 
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Figure 1: The closed timepath contour C + © C_ 



Ptf [<f> f ] = J v^ptM J p<M¥- 8[</>+{t f ) - $/] ex P {^s[0 + ] - sy>-]) 

(17) 

where 5[0 + (i) — $/] is a delta functional, imposing the constraint + (x, t) = $/(x) 
for each x. This is no more than the statement that, for a given initial state, 
the probability amplitude is given by the integration along C + , and its complex 
conjugate (which, when multiplied with it, gives the probability) is given by the 
integration back along C_. The ± two-field notation is misleading in that it suggests 
that the + (= a ,+) and 0_ fields are decoupled. That this is not so follows 
immediately from the fact that 0+(t/) = 0_(i/). 

To return to the initial conditions, the simple analytic results of the previous 
section could only be achieved if Pt[&] were Gaussian and, therefore, that p t() [$] 
itself be Gaussian. This may not be unrealistic. The simplest such distribution has 
$ Boltzmann distributed at time to at a temperature of To = Pq 1 and chemical 
potential fi according to a /ree-field Hamiltonian Hq, which we take as 



H n = / d 3 x 



18) 



where ir a = a in real time. We stress that Gaussian does not necessarily mean 
free. Most simply, particles need to interact before they can equilibriate. More 
importantly, the free-field Gaussian approximation adopted here can be extended 
to include interactions self-consistently in a Hartree approximationS. 

Since chemical potentials are not usually relevant to relativistic bosons a few 
words are in order. The 0(2) invariance of H Q leads to a conserved Noether current 
with conserved charge 

Q = J d 3 X (02^1 - 7T20l) 

d 3 x (0201-0201). (19) 



The numerical value of Q is the number of particles minus the number of antipar- 
ticles. The thermal probability distribution p to [$] is thus taken to be 

p t M = (^^ \e~^ H °-^,t ). (20) 

We note that interactions with weak coupling will not change pt [$] significantly. 

There are two ways to proceed. The first, which to us is the most natural, 
accepts H as determining the temporal evolution of the fields, and relegates fi to 
the initial boundary conditions on the fields. From this viewpoint, Pt [$] can be 
written as the imaginary-time path integral 

PtM= [ 2>0exp(zS o [0]), (21) 
for a corresponding action 



SM] = d 4 x 



^(^0 a )(^0 a )-im 2 02 



(22) 



The boundary condition -B M [$] incorporates the chemical potential. In terms of the 
Q eigenstates, the complex field — (0i + z02)/v2 an d its adjoint, -B M [$] is 

B^} : 0(t o ) = $ = e-^0(*o - ift), (23) 

where q — 1 is the 0-field eigenvalue of Q and the time-integral in Eq.|22] is taken in 
imaginary time from i to to ~~ ifla- I n terms of the same complex fields we have 



S Q [<t>] = J d 3 x[(<%0*)(<9» - m 2 0>] (24) 

and, with tt = (tti — iff^/y^, 

H = J d 3 x[n*7T + (V0*)(V0) + m 2 0*0]. (25) 

However, for calculating vortices a decomposition in terms of <f) a , n a is usually prefer- 
able. 

Although we are just setting an initial condition the effect is, inevitably, to give 
an action So[0] of the form of S[<j)} of Eq.[l6|. This permits the interpretation that 
the action S[<p] is valid for all times t, with the proviso that the system is in thermal 
equilibrium for t < to, during which period the mass m(t) takes the constant value 
m and \(t) = 0. 

On relabelling the integration variable of Eq.^T]by 03, we now have the explicit 
form for p tf [$/]:- 

p [$-] = f V® [ £>0 3 e i5o[<fel [ V(j) + <$>V(j)_e l{s[(t,+] - s[4, - ]) 5[<i) + {t f ) - $,] 

VfaVt+Vt- exJiSolfa] + i(S[<f>+] - S[<l>-])\6[<f> + (tf) - $/], (26) 



B, 



where the boundary condition is now (in terms of the field combinations — 

(0x + z0 2 )/V2) 

B, : Mto) = e-^ q Mto - ifo). (27) 
More succinctly, p tf [$ /] can be written as the time ordering of a single field doublet:- 

p tf [$ f ] = I V(j>e iS ^ 5 [(f> + (t f ) - $ f ], (28) 

J Bp 

along the contour C = C+@C- ©C3 of Fig. 2, extended to include a third imaginary 
leg, where takes the values <f> + , 0_ and 03 on C+, C_ and C3 respectively, for 
which Sc is 5'[0 + ], S[(j)-] and £[03], for which last case m(t) = m, \(t) = 0, t G C3. 

As a final step in these formal manipulations we see that expression Eq.|2^ en- 
ables us to write the $-field ensemble averages (...) t in terms of the 0-field Wightman 
functions. Consider the generating functional:- 

Zp\j+,j-,j 3 ] = J b V(j> exp jzS c [0] + i J j0|, (29) 



to 
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Figure 2: A third imaginary leg 



where / jcf) is shorthand for:- 

j<j> = J o dt [j + (t)<P + (t) ~ + J Q h(t)Mt)dt, (30) 

omitting spatial arguments. Then introducing a a (x) where a = 1,2, we find:- 
p t/ [$] = J Va jf V<}> exp|z5 c [0]| expjz J cPxai a (ic)[<P + (t f , x) - $(x)] a } 

= fva expj-z | a a $ Q |z M [a,0,0], (31) 

where a is the source a?(x, £) = a(x)5(t — £/). As with £>0, Va denotes Y\i Va a . 
Ensemble averages are now expressible in terms of Z^. Of particular relevance, 

W ab {\v-v'\-t) = ($ a (r)$ 6 (r')) t 

= (0 a (r,t)0 6 (r',t)}, (32) 

the equal-time Wightman function with the given thermal boundary conditions. 
Because of the time evolution there is no time translation invariance in the double 
time label. 

Not surprisingly, pt[$] can only be calculated explicitly in very simple circum- 
stances, like our Gaussian approximation, but before we do that we still have to 
extract the relativistic limit. Since the chemical potential is embedded in the equi- 
librium boundary conditions and, like temperature, can only be defined in equilib- 
rium, this is essentially an equilibrium problem and, for the moment, we forget the 
dynamics. 

3.2. Deriving the N on- Relativistic Limit 

The expression Eq.^8| is valid for all \i and all To but, as it stands, is not sympa- 
thetic to the isolation of a non-relativistic limit in a form with which we are familiar. 



Nonetheless, the extraction of a non-relativistic regime from it is not difficult. In- 
stead of working with p to [$] directly, we can work with the partition function Z^, as 
we saw in the previous section. The partition function for this equilibrium theory 
is (for doublet sources j a restricted to the CVcontour) ^[0,0, j'3] of Eq.|29|, written 

Z^j] = Z M [0, 0,j] = j B V<\> exp{iS Q [0] + % J j a a }, (33) 

where «So[0] is given in Eq.^2| and we integrate only along the contour C3. As before, 
"D0 ee Z>0]Z>0 2 - On rotating to Euclidean time r = it, Z^j\ becomes 

ZM = (exp{| ja<Pa}) = J b V<\> exp|-5o, s [0] - J J a a }, (34) 



where Sq eI^} is the relativistic Euclidean action 



Sq,e[4>\ = j dr I d 6 x 



(35) 



and the sum is taken over fields = (0i + 202) / v2 satisfying the boundary condition 
B^: 0(x, r) = e~^ oM 0(x, r— /? ) i n imaginary time. The dot now means differentiation 
with respect to r, a = <9 r a . 

The bracket (...) here denotes the thermal average 

(F[(j)}) = Tr{e~ MHo ->* Q) F[0]}. (36) 

where Hq is given in Eq.[18|. Since it is time-independent we have dropped the 
t-suffix. 

Zfj\j\ can be calculated as it stands. However, because of the chemical potential 
in the boundary conditions it is not immediately obvious how to identify the phase of 
the system. This is clarified by adopting a second approach to chemical potentials, 
in which we transfer the chemical potential term fiQ to the Hamiltonian, to create 
an effective Hamiltonian H — fiQ. Although H — \iQ no longer generates time 
translations of the physical fields and 0*, it does generate time translations of the 
effective field <fi(t) = e tfJ,t 4>(t) (and its adjoint). To see the advantages of removing 
\x from the boundary conditions we express Z^ of Eq.|33| not as a sum over 0-field 
histories but as a sum over 0-field histories. Although not usually posed in this way 
the details are well-jmderstood. For example, see Kapusta0, Haber & WeldorMI, 
and Bernstein et alS 

With the Jacobian unity, Z^j] becomes 

Z,[j] = J P0 exp{-S M [0] - J jj a y (37) 
where S M [0] is now the effective action 



rPo , ., 
S^} = J dr I f /•'./■ 



\4>l + ^(V0 a ) 2 + \ml4>l + ^(0201 - 2 0l) 



(3? 



in which 



2 2 2 

m = m — n 



(39) 



It follows from Eq.[23| that, with \iQ bundled with H the <p a integrations are now 
taken over periodic configurations with boundary conditions B Q : a (x, r) = a (x, r— 
f3o) in imaginary time, with no /i-dependence. Most importantly, the zeroes of <f) 
are the zeroes of 0, leading to identical line densities. 

This displacement of \i from the boundary conditions to the action has enabled 
us to replace the classical potential by an effective potential in which the state of 
the system is more transparent. Semiclassically, when rriQ > (i.e. fi 2 < m 2 ) the 
0(2) symmetry is unbroken. This describes our disordered initial state. However, 
once j«q < (i.e. /i 2 > m 2 ) the free theory is unstable. This is a signal that the 
0(2) symmetry is broken and a transition to an ordered phase has occurred. Thus, 
effectively, it is m 2 ,, rather than m 2 that carries the time dependence, changing from 
positive to negative at t = t . 

Semiclassically, the non-relativistic limit is trivial. Although we achieve it by 
following the phase boundary m 2 , = from high T and low fi to low T and high 
H the decomposition of m. 2 , is immaterial. Only in the cross-terms in(<p24>i — 4>2 ( j ) i) 
will the change be noticed. Even this is misleading since, as we shall see, this term 
plays no role in our Gaussian approximation. 

All that remains is a problem of interpretation. Let us first consider the rela- 
tivistic regime in which the initial environment is very hot, with T ^> m ^> fi. In 
this case, with fi irrelevant, the symmetry is broken by a change in m 2 . Suppose 
that, on the completion of the transition, m 2 , ~ m 2 takes the value m 2 , = — M 2 < 0. 
If, in this relativistic case, the final temperature is very low, there are then no 
thermal effects once interactions have been taken into account and M is a physical 
parameter, determining the Higgs mass, m H = \/2M. However, because of thermal 
radiative effects m (t < 0) = mo is not a physical parameter. In practice, there is no 
loss in taking it as the effective scalar field mass at temperature T = (3q 1 , despite 
the fact that such a mass is defined in terms of large scale fluctuation averages (the 
effective potential). That is, in the mean-field approximation, 

m 2 = -M 2 (l-^). (40) 

The results will turn out to be largely independent of mo provided it is comparable 
to M, which is the case if we do not quench from too close to the transition. On the 
other hand, for the non-relativistic theory, for which T/m 1, the parameter m can 
always be identified with the boson mass. In this case we change the sign of m 2 , in 
the action by a change in fi from /i 2 < m 2 to fi 2 > m 2 , a density quench. To preserve 
uniformity of notation, we take m^(t) = uIq when t <t and m^t) = —M 2 < once 
the transition is complete. 



3.3. Extracting the Non-Relativistic Field Theory 



At this stage we could proceed directly to our Gaussian fluctuation model. How- 
ever, it is advantageous to show how the unfamiliar equation Eq.|3^ can be cast in a 
form standard to non-relativistic many-body theory. In fact, one of the difficulties of 
making straightforward comparisons between early universe physics and many-body 
physics is this difference in formalism, a reflection of the freezing out of antiparti- 
cles in the non-relativistic regime. Thus, for example, a non-relativistic order field, 
with its interpretation as a probability amplitude for finding particles, is necessar- 
ily complex. As is apparent, our main tools in relativistic particle physics are path 
integrals which, conventionally, are not separated into particle and antiparticle sec- 
tors. This split has to be made if we are to recover the familiar non-relativistic as 
a field theory in its own right. 

To do this we rewrite Z M of Eq.^7| as a sum over histories in the full phase-space 
of the fields (p a and their conjugates 7f . As we shall see, it is just this extended sum 
over histories that is needed for isolating the non-relativistic limit of the theory. 

As a phase-space integral, Z^[j] takes the form 

ZJuX = J b V^V^V^V^exp^-S^n,^} - J jj a }, (41) 



where 



-S^ir, 0] = r° dr{ ( / d 6 x m a <p a )-[H --fiQ)\. (42) 



o 



We have made the integration over the four field degrees of freedom explicit. Al- 
though the <f> — a-fields are constrained to be periodic, the n a integrations are un- 
constrained, and S^e cannot be identified directly with the 'effective' action in the 
presence of a chemical potential. However, the tt integrations are readily performed 
to recreate of Eq. [37[ 

It has been shown in the work of Kapusta, Bernstein et al.@ how, in the limit 
m ^> T ~ — /i nr , where fi nr = fi — m, we recover the energy spectrum (and free 
energy) of a non-relativistic Bose gas. We shall now show that this same phase- 
space analysis permits us to recover the familiar non-relativistic partition function, 
for free bosons, 



bo 



Z,JV*,V] = / VrVi, exp -S,jr^] + / (#* + , (43) 



where the boundary condition 6 denotes simple periodicity of the fields in 
imaginary time with period (3q. In Eq.|l3| S ^[ip* , ip] is the non-relativistic action 



* ' 3, 



kri - r^p) + ^(w>*)(w) - fi m r^ 

2 2m 



o 



dr I d 3 xil> 



— - —V 2 - 

dr 2m ^ T 



if) (44) 



in the order field ip and its conjugate. 

The non-relativistic field ip(pc,t) is understood as annihilating a particle at the 
point x, whereas ip* (x, t) creates a particle at x. Although the particles are distinct 
from the antiparticles there are no antiparticles in Z^. \ip\ 2 measures the local 
particle density, a role not available to |0| 2 of the relativistic field theory. Of course, 
we can decompose ip as ip = ^75(^1 + ^2) but, for once, we decline to do so since 
ipi and ip2 have no obvious interpretation. 

Zfj\j\ of Eq.|37] and Z^ ni [rj*, 77] of Eq.^3| cannot be identical, coming from different 
origins, but they are closely related in the regime m ^> T ~ — fi nr . To see this, we 
rewrite Z^j] of Eq.^T| in terms of the complex fields 0, <f)*, it, n* as 

Z^f] = J P7r*X?7rP0*P0exp{-^[7r,0] - J (f(f> + j<f>*)}, (45) 
where we have dropped all tilde labels from the fields and 

- -S^e[kA] = J^ dT {J d " x + - (#0 - -V0)\, (46) 

in which H is given in Eq.[25| and Q is now, equivalently 

Q = i f d 3 X (0*7T* -7T0). (47) 



We could proceed as before, and perform the tt, it* integrations, which would 
give us Zplj] of Eq|J7] in complex field form, but this is nothing new. Instead, we 
perform canonical transformations from the phase-space fields 0, 0*, ir, tt* to a new 
set of fields that, in the non-relativistic limit, annihilate and create 

particles and antiparticles. Specifically, 

** = i/r^-vsr- (48) 

annihilate and create non-relativistic particles in this limit, whereas 



2 V2m 

= &_ * tt*, (49) 

annihilate and create non-relativistic antiparticles respectively. 

To confirm that this is the case we invert the above. The field measure 



(50) 



is unchanged. Inserting these decompositions in Z^lj, j*} of Eq.|35| (where, for sim- 
plicity, we switch off the sources j,j*) gives of the form 

Z„ = ^[0,0]= / VWVVV&VVexpi-S^E} (51) 

J Bq 

where S^e of Eq.^6| is now a functional of etc., but i?o is a boundary condition 
on the cf) fields, and not the \I> fields. 
Then, on direct substitution is 



o 



— - —V 2 + (m - /i) 

or zm 



+ jj dr j d 3 x J- j(Vtf)(W) + (W*)(W*)j, (52) 

after removing total derivatives. This can be written in terms of the non-relativistic 
action S^ip*, t/>] of Eq.g| as 

S,,e = <W** , *] + S_ M _ m [**, *] + So./t**, * , **, (53) 

where Soj, coupling particles to antiparticles, is the third term of equation Eq.|52| 
and n nr = fi — m. S^^Sf*, \P] and S'-_ A1 _ m [ 1 l , *, ^] each possess an 0(2) invariance, for 
the individual conservation of particle and antiparticle number. The 0(2) symmetry 
of the antiparticle action S-n- m [^*, ^f] remains unbroken. It is the 0(2) symmetry 
of the particle action of Eq.[|4] which leads to non-relativistic vortices on its 
breaking. Away from the non-relativistic limit the 0(2) ® 0(2) symmetry of S^e 
is reduced to the 0(2) of the relativistic theory (the conservation of particle minus 
antiparticle number) by So,/- 

To see this pattern we observe that, on partial integration, Sqj can be thought of 
as providing sources for ^ and ^f* in of Eq.^T], which can be integrated out. This 
gives Z^ as an integral over the non-relativistic partition with chemical potential 
— m — (jl, 

Z^= f D$*P$ expi-S.J**, $]} Z_ m _ /1 [-V 2 */2m, - - V 2 **/2m] (54) 

J Bo 

in the terminology of Eq.|i3| where, again, the periodic boundary condition is on 
the ip fields. For the free-field case in hand, Z_ m _^ in the integrand is Gaussian 
in its sources. For the non-relativistic regime, in which — m\ <C m, whence 
—fi—m « —2m, the factor Z_ m _ M « Z_ 2m in the integrand of Eq.[53]can be expanded 
in powers of momenta. For particle kinetic energies k 2 /2m <C m, Z_ rn _^ « 1. For 
typical k = |k|, with k 2 /2m = 0(T), the errors are 0(k 2 /m 2 ) = 0(T/m). From 
Eq.|4£| and Eq.|4^, we identify the semi-classical with \I/ as 

= + ~ - 1 =^> (55) 

V2m v2m v2m 



on neglecting the anti-particle field fluctuations. As a result, zeroes in <fi are also 
zeroes in $. Further, on neglecting these fluctuations, <fi ~ n* /2m showing that, in 
the non-relativistic limit, tc is approximately periodic in imaginary time. Thus the 
boundary condition B , denoting periodic fields, becomes the boundary condition 
bo for periodic \I> fields, whence we achieve our goal of showing that, in the non- 
relativistic limit, the relativistic partition function 

Z, « = [ P**P*exp{-^ nr [**, vl/]} , (56) 

where fi nT = fi — m as before. The inclusion of sources will not change this result, 
although we shall have to include sources that couple to the tt*, ti fields. 

This demonstration of how calculations with the relativistic of Eqj45| are 
indistinguishable from those of the non-relativistic when /i »s m can be extended 
to the case of symmetry breaking, although we shall not be able to go beyond the 
simplest Gaussian approximation in this talk. Further details are given in our 
workEJ. 

Now that we have demonstrated how to recover the non-relativistic field theory, 
we return to the relativistic formalism of Eq.|3S], for large and small \x. 

4. A Gaussian Model for Vortex Formation 

Having established the role of the chemical potential in letting us interpolate 
between relativistic and non-relativistic regimes in the initial conditions (with a 
corresponding interpolation between relativistic and non-relativistic field theories) , 
we are now in a position to determine the effect of these initial conditions on vortex 
production in a simple model of Gaussian field fluctuations. 

We have already assumed that the initial conditions correspond to a disordered 
state in equilibrium for t < to- The action Eq.|T^ was proposed in a context in which 
the chemical potential was relegated to the boundary condition. As we observed 
earlier, if we wish to represent it so as to match the action Eq.|38], in which the 
chemical potential is manifest, we can accomodate changes in m or changes in fi 
by the requirement that, at t = t , rn^{t) = m 2 {t) — fi 2 (t) changes sign. Most 
simply, we assume that, for t > to, m^i) takes the negative value m^(t) = —M 2 < 
immediately. That is, the potential at the origin has been instantaneously inverted 
everywhere, breaking the global 0(2) symmetry. If \(t) = A is very weak then, for 
times Mt < ln(l/A), the 0-field, falling down the hill away from the metastable 
vacuum, will not yet have experienced the upturn of the potential, before the point 
of inflection. Thus, for these small times, X(t) can also be set to zero, and is 
Gaussian, as required. 

In the relativistic regime /i -C m, and the t-dependence of m is entirely carried 
by m, as in Eq.[iO|. On the other hand, in the non-relativistic limit, when /i ~ m 
and m is unchangine in time, wig changes from 



ml(t) = )71q = —2mfi > 0, 



(57) 



say, when t < to, to 

m 2 Q (t) = —M 2 = -2mfi f < 0, (58) 

when t > to- Henceforth, we take to = 0. 

The onset of the phase transition at time t = to is characterised by the in- 
stabilities of long wavelength fluctuations permitting the growth of correlations. 
Although the initial value of (0) over any volume is zero, we anticipate that the 
resulting evolution will lead to domains of constant (0) phase, whose boundaries 
will trap vortices. 

Of course, instantaneous change is physically impossible. Consider small ampli- 
tude fluctuations of a , at the top of the parabolic potential hill. Long wavelength 
fluctuations, for which |k| 2 < M 2 , begin to grow exponentially. If their growth 

rate Q(k) = \J M 2 — |k| 2 is much slower than the rate of change of the environment 
which is causing the quench, then those long wavelength modes are unable to track 
the quench. It will turn out that the time-scale at which domains appear in this 
instantaneous quench is t^ = 0(M~ 1 ). As long as the time taken to implement the 
quench is comparable to t^ and much less than tf = 0(M~ 1 ln(l/ X) the approxima- 
tion is relevant. We note that, in the non-relativistic regime, for which Q(k) has 
the same definition, 

n 2 (k) = M 2 - |k| 2 = 2m(ii } - ^) 

= 2m(fx f - e(Jfe)). (59) 

Thus the momentum restriction |k| < M is just e(k) < [if. 

We are now in a position to evaluate pt[$\, or rather W a b(r;t), for t > 0, and 
calculate the defect density accordingly. Before we quote the result we note that the 
i/i(0 2 0i — 2 0i) term in S^cp] of Eq.|38"| couples the a = 1 and a = 2 fields (fi a together 
and, in general, G a fe(r — r'; t, t') = (0 o (r, t)cj) a (r', t')) is not diagonal in the 0(2) labels. 
However, for equal times diagonal behaviour is restored as G a b(r;t,t) = 5 a t,G(r;t,t) 
and the i//(0 2 0i — 020i) term can effectively be discarded. This leaves us in a 
situation for which Eq.|^ and Eq.[Tt]| are satisfied and the results of Halperin are 
directly applicable. That is, W a b is diagonal, 

W ab (\r -r'\;t,t)= 5 ab W(\r -r'\;t,t), (60) 

whence W(\r — r'\;t) = (<p a (r, t)(p a (r', t)) (no summation), the thermal Wightman 
function for either (pi or 2 . 

But for the chemical potential, this situation of inverted harmonic oscillators 
was studied many years ago by Guth and PiliSl and Weinberg and WiiH In the 
context of domain formation, we refer to the recent work of Boyanovsky et alHI, 
and our owni'0. For the case in hand, if we make a separation into the unstable long 
wavelength modes, for which |k| < M, and the short wavelength modes |k| > M, 
then W(r; t) is 

W(r;t) = f fke ik - x C{k-fi)\l + A{k){cosh(2fl{k)t)-l) 

J\k\<M I 



+ 



|k|>M 



1 + a(k)(coa(2w(k)t) - I) 



(61) 



with r = Ixl and 



n 2 (k) 

w 2 (k) 
A(k) 

a(k) 



M 2 - Ikl 2 , 



-M 2 + Ikl 



1 + 



in 



uo 2 {k) 

n 2 (k) 

u 2 {k) 
w 2 (k)J 



(62) 



All the /i-dependence is contained in the factor 



2u{k) 



where 



uj 2 (k) = |k| 2 + m 2 . 



(63) 



(64) 



In the zero chemical potential limit /i = m + /io = 0, C(k,fi Q ) takes the familiar 
form 

C(fc, /io) = ^yC0t/i(^(A;)/2) (65) 

On the other hand, in the non-relativistic limit e(k) = k 2 /2m ~ u m and /io <C m 

C(k,no) is the equally familiar Bose distribution 



C(fc,/x ) 



2u;(fc) VI - e -^( e ( fe )-w) /' 



(66) 



Further, in the high temperature relativistic limit T 3> m Q , C(k,fi ) of Eq.^5] sim- 
plifies as 

= ^ (67) 

Equally, in the non-relativistic limit when e(k) — /i <C T , 

To 



oc 



e(fc) + |/i | 

To 
k 2 + mo 



(6* 



on using the definition of m given earler in Eq.|58|. From Eq.|13] onwards it follows 
that the overall scale of W is immaterial to the vortex density. As a result, in these 
regimes the relativistic and non-relativistic W(r,t) are identical, once we take the 
identifications of Eq.[58] into account. 



We observe that, if we were to use W(r; t) of Eq.|61] as it stands, then both 
W(0;t) and W"(0;t) necessarily suffer from ultraviolet divergences. However, the 
string thickness at the end of the quench will be 0(M~ 1 ). It is the zeroes coarse- 
grained to this scale that will provide the subsequent network. Thus, if we do not 
probe the field zeroes within a string we need consider only the first term in Eq.[5T]. 
There is another point. Even before the quench begins there is a high density of line 
zeroes coarse-grained to this same scale 0(M~ l ) in the initial equilibrium phase. 
However, these modes are entirely transient. If we were to calculate the correlations 
of Pi(x) at differenttim.es t and t', we would find rapidly oscillating behaviuor with 
period At = 0(m -1 ). On the other hand, a calculaation of the density correlations 
at different times from the unstable modes in Eq.^l] does not give oscillatory, but 
damped, behaviour. It is the residue of the strings produced by the unstable modes 
that survives to produce the network, and the transient strings can be ignored. 
Henceforth, we retain only the first term 



W(r;t) = / fke lk - x C(k,fM)) 

J\k\<M 



1 + A(fc)(cosh(2fi(fc)t) 



(69) 



of Eq.^T[ In the two critical regimes where Eq.^JT] and Eq.|6^ are valid, W(r; t) is 
again the same for both cases. 

Even though the approximation is only valid for small times, there is a regime 
Mt > 1, for small couplings, in which t is large enough for cosh(2Mt) w |exp(2Mt) 
and yet Mt is still smaller than the time 0(lnl/A) at which the fluctuations begin 
to sample the ground-state manifold. In this regime 

W(r;t)~ [ fkC{k^y^A(k)e m{k)t . (70) 

J\V\<M 

In these circumstances the integral at time t is dominated by a peak in the integrand 
fc2 e 2n(k)t a ^ arounc i w here 

** = M ( 1+0 (-m))' (71) 

and we have assumed M and tuq to be comparable. The effect of changing (3q is 
only visible in the 0(1/Mt) term. In fact, we are being unnecessarily restrictive in 
wanting to preserve the identical behaviour of Eq.^ and Eq.^Sl All that is required 
for identical leading behaviour in relativistic and non-relativistic regimes is that 
C(k, /io) varies slowly in the vicinity of the peak of the integrand at tk 2 c ~ M. [A(k) 
is already slowly varying]. This is the case when, allowing for coefficients 0(1), 



K ^ = K«tM<ln(\) (72) 



where, in the same spirit, we have taken |/xq| < yU/- The upper bound on tM is 
a reminder that interactions are always present, and the Gaussian approximation 



must fail as soon as the field fluctuations have extended to the true ground-states at 
the minima of the potential. The lower bound is necessary for the integrand to be 
peaked strongly so that the saddle-point approximation is valid. As long as Eq.[7^ is 
basically correct, any difference between the relativistic and non-relativistic regimes 
will be non-leading. 

Assuming these limits, we recover what would have been our first naive guess 
for a correlation function based on the growth of the unstable modes, 

W(r;t) ~ / 4 3 ke^ r e m{k) \ (73) 

J\k\<M 

We understand the dominance of wavevectors at k 2 = M/t in the integrand as 
indicating the formation of domains of mean size 



^(t) = 0(y/i/M) (74) 



once Mt>l. Specifically, we take f (t) = 2^/t/M. Once Mt > 1 then £(t) > M~ l , 
where M -1 characterises the cold vortex radius. In the weak coupling approximation 
individual domains become large enough to accomodate many vortices before the 
approximation breaks down. There is no difficulty with causality since domains 
increase in size as £ = ^7= < 1. On neglecting terms exponentially small in Mt, 
W(r; t) of Eq.[73] can be further rewritten as 

roc „ , 

W(r; t) ~ e 2Mt / dk sinc(Ar) k 2 e~ tk /M (75) 
Jo 



rp2 



W(0;t)exp[-— -) (76) 



where 



e 2Mt 



w{a ' t] " c WW (77) 

for some C. The exponential growth of W(0; t) in t reflects the way the field am- 
plitudes fall off the hill centred at $ = 0. With the peaking in wavelength I = k~ x 
understood as indicating the appearence of domains of characteristic linear dimen- 
sion the Gaussian in r is a reflection of the rms variation A£ in domain size £. 
This variation is large. If we isolate the Gaussian saddle-point in Eq.[73| as 



W(r; t) ~ e 2Mt / dk smc(kr) kie~^> /2 ^ k > (78) 
Jo 

then 

^ = ^ = i (79) 

To calculate the number density of vortices at early times we insert the expres- 
sion Eq.|75| for W into the equations derived earlier, to find 

..11 



We note that the dependence on time t of both the density and density correlations 
is only through the correlation length £(t). We have a scaling solution in which, 
as the domains of coherent field form and expand, the interstring distance grows 
accordingly. Since the only way the defect density can decrease without the back- 
ground space-time expanding is by defect-antidefect annihilation, we deduce that 
the coalescence of domains proceeds by the annihilation of small loops of string. 
However, because the density of vortices only depends on £(t) in this early stage, 
the fraction of string in 'infinite' string remains constant. Thus, at the same time 
as small loops disappear, other loops must rearrange themselves so that the length 
of 'infinite' string decreases accordingly. Finally, there is roughly one string zero 
per coherence long held belief for whatever mechanism. 

There is one final concern. A necessary condition for this rolling down of the 
field to be valid is that the initial field fluctuations should be small enough that 
there is no significant probability that the field is already in the true vacuum. As 
long as we do not begin the quench from too close to the transition, there is, in 
fact, no difficulty. See our work0 for more details. 



4-1. Vortex Density Correlations 



In addition to the gross vortex density Eq.|80| we can calculate the density-density 
correlation functions CV,(r; t) of Eq.[l^, identical for both the relativistic plasma and 
the non-relativistic medium. 

Yet again, as in the case of the density n{t), the t-dependence of Cy only occurs 
implicitly through £(;£). The simple analytic form of Eq.[F5] enables us to calculate 
A and B, up to exponentially small terms in Mt, as 



A(r;t) 
B(r;t) 



-2r 2 /e(t) 



n 2 ^(t) (l-e-^W)) 2 



-2r 2 /£ 2 {ty 



~2 n 



< o, 



7T 2 e 4 (t) (l-e-2- 2 /C 2 W) 
Suppose that r = (0, 0, r). Then 



> 0, 



*1) 



C^^ — C\ 



22 



A(r;t) - B(r;t) 



2r 2 



,-2r 2 /? 2 (t) 



n 2 £ 4 {t) \£ 2 {t)J {l-e- 2r 2 /* 2 (t))s 



< 0. 



52) 



That is, we have anticorrelation of densities for parallel directions (and positive B 
for orthogonal ones). 

It is useful to expand Eq.|81| and Eq.|83 for r < £. On normalising by a factor of 



n- 



A{r;t) 
n 2 (t) 



-l + O 



< 0, 



B(r;t) f e(t) 



n 2 {t) \ r 



,2 



0(1) > 0. (83) 



For r > £ there is exponential falloff but, from Eq.|83| we see that, in units of n(t)~ 2 , 
the anticorrelation is large. Since strings with a long persistence length would imply 
positive parallel correlations, we can interpret these anticorrelations (and others) as 
a reflection of an increased string bendiness. Although it is difficult to be precise, 
this suggests a significant amount of string in small loops. This is an important 
issue since, as we noted, early universe cosmology requires infinite string if string is 
to be the source of large-scale structure. In practice, some is enough. 

There is another, indirect, way in which we can see the tendency for more string 
to be in small loops than we might have thought. In numerical simulations of string 
networks, the rule of thumb for a regular domain structure in field phase is that a 
significant fraction of string, if not most, is infinite string. 

However, as we saw earlier in Eq.[79|, we do not have a regular domain structure 
in our model but have domains with a large variance A£/£ = |. Unfortunately, the 
domain structure that we have here is not yet appropriate for a direct comparison 
since, as well as the domain size, the field magnitude has a variance about = 
0{Me Mt ). The work of Guth and Pi@ shows that this can be parametrised by an 
effective dispersion in t in of At — 0(M~ 1 ). Nonetheless, consider an idealised 
case in which domain growth stops instantaneously because of back-reaction at some 
time tf. The distribution of strings will then be as above for £ = £(t/), while the field 
adjusts to the vacuum manifold, without changing phase, in each domain. We would 
expect some string-antistring annihilation to continue while this adjustment occurs, 
so that the n{t) calculated previously is an overestimate of the string density at the 
end of the transition. However, the domains will still be of varying size with variance 
still given approximately by Eq.[79]. The inclusion of domain variance in numerical 
simulation of string networks (albeit in a different way from that proposed here) 
shows" that, the greater the variance, the more string is in small loops. Beyond 
observing that there seems to be some infinite string, we will say no more. 

All our results are for quenches that go from significantly above the transition to 
significantly below it. These give the smallest value of A£/£ possible and, plausibly, 
the best chance of producing infinite string. Had we begun closer to the transi- 
tion, where the initial fluctuations are larger, the dispersion in domain size, and 
consequently the fraction of string in small loops, would be even greater. 



5. Conclusions 



In this paper we have shown how global 0(2) vortices appear, at a quench from 
the ordered to disordered state, as a consequence of the growth of unstable Gaussian 
long wavelength fluctuations. Most importantly, in the light of discussions about 
the extent to which vortex production in low-temperature many-body systems sim- 
ulates vortex production in the early universe, our model supports the analogy. 
Specifically, with our simple assumptions, vortex production is identical in both a 



relativistic high-temperature quench, in which the initial state is characterised by 
T ^> m, and in a non-relativistic density quench in which the initial state is de- 
scribed by T C m (to freeze out antiparticles) with a chemical potential /i nr ~ T. 
All that is required is an appropriate translation of the parameters from the one 
case to the other. To see this, we have shown how to isolate the non-relativistic 
particle/antiparticle sectors of the relativistic path-integrals that we use in our cal- 
culations. The 0(2) symmetry, whose breaking leads to the production of vortices, 
leads to the conservation of particle number minus antiparticle number in the former 
case, and to the conservation of particle number in the latter. 

In our simple model of Gaussian fluctuations the resulting string configura- 
tions scale as a function of the correlation length £(t) = 0(t?), at about one vor- 
tex/correlation area. This is compatible with the Kibble mechanism for vortex 
production on domain boundaries upon phase separation. For a weak coupling the- 
ory the domain cross-sections are significantly larger than a vortex cross-section at 
the largest times for which the approximations are valid. However, there is a large 
variance in their size, with A£/£ = |. Because of this there is more string in small 
loops than we might have anticipated. 

We stress that our model can, at best, describe weak coupling systems for the 
short times while the domains are growing before the defects freeze out. This is 
unsatisfactory for most early universe applications and for low-temperature many- 
body systems. However, we know in principle^ how to include back-reaction (still 
within the context of a Gaussian approximation) to slow down domain growth as 
the field fluctuations spread to the ground-state manifold. The identity of the 
weak-coupling results of the two regimes should survive to this case also, although 
it will probably lead to different conclusions from those above. Further, there is no 
difficulty in principle of embedding these results in a FRW metric, along the lines 
oie-a. This is being actively pursued. 
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